POWER SYSTEM STATE 

ESTIMATION : IMPLEMENTATION AND EVALUATION 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

P. M JALTARE 


to the 

DEPARTMENT OF ELECTRICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

AUGUST, 1978 



Dedicated to 


AAJOBA 



! ’tri^ 

-7- u ijini. 

4 mw Nob Jnk 

’ 1 AUG 1978 


~fL- 

7 5\g f 


& 9 -~ ^ ~ 



OERTIPICATE 


Certified that the work entitled "Power System 
State Estimation: Implementation and Evaluation", hy 
Mr. P.M. Jaltare has been carried out under my 
supervision and has not been submitted elsewhere for 
a degree. 


August 1978. 






Pr. 


R. P , Aggarwal 
Professor 

Department of Electrical Engineering 
Indian Institute of Technology 
Kanpur 



AOKKOWLEDGMENT 


I express my sincere tharks to Er. R.P. Aggarwal 
whose valuable guidance and cooperation has encouraged me 
at every step. I feel very much obliged to Dr.M.A. Pai 
for his keen interest and continued help in this topic. 

I also thank all my friends, '0 Block' friends in 
particular who helped me directly or indirectly in carrying 
out this work and made my stay here more enjoyable. I thank 
Vishwas Eeshpande for the help he has been doing to me from 
the beginning of our stay here. Thanks are also due to 
Mr. K.IT. Tewari for the fast typing he has done. 


P. M. Jaltare 



iv 


OONTETOS 

Page 

CHAPTER 1: INTRODUCTION 1 

1.1 Motivation 1 

1.2 Summary of "Various Chapters 

3 

chapter 2: LTOIPED ¥LS AND EAST DECOUPLED STATE 

ESTIMATION 5 

2.1 Statement of the State Estimation Pro'blem 5 

2.2 The Lumped WLS Estimator 7 

2.3 Model for Pover System Network 9 

2.4 Entries of the Jaco"bian 12 

2.5 Implementation 15 

2.6 The East decoupled estimator (EDS) 18 

2.7 Implementation 20 

2.8 Acceleration of Convergence 22 

2.9 State Estimation Based on Real Power 
Plow in LinevS and Voltage Magnitudes 

only (PV) 23 

2.10 Implementation 24 

CHAPTER 3: THE LINE ONLY ALCOUTHM (»L0' aLC0RITH14) 26 

3.1 The AEP Line Only Algorithm 26 

3.2 Modification of AEP Algorithm to take 

into account efficiently both End Line 
Plows 29 

3.3 Implementation 31 

3.4 Computer Results 34 



CHAPTER 4: STATE ESTIMATION BY INI" -PENDENT 



EQUATION SOLUTION 

35 

4,1 

Splitting lumped Equations into 
Independent Equations 

' 35 

4.2 

Solution Procedure 

38 

4.3 

Solution for Voltage from Line Plows 

40 

4.4 

Computer Results 

42 

CHAPTER 5: 

DETECTION AND IDENTIPICaTION OP CROSS 
MEASURMENT ERRORS 

44 

5.1 

Statement of the Problem 


5.2 

Mathematical Baclcground 

45 

5.3 

Implementation 

50 

CHAPTER 6: 

DYNAMIC STATE ESTIMATION 

53 

6.1 

Application of Observer Theory to 

Pov/er Systems 

55 

6.2 

Extended Linear Observer 

55 

6.3 

Implementation 

58 

CHAPTER 7: . 

“’'UlTHER^RES^i.RCH-TOPrCS IN STaTE 

59 


ESTIMATION 

» 

^ 7.1 

Topics Related to State Estimation 
Applications 

59 

7,2 

Topics Related to Improvement in 
Efficiency of Static State Estimators 

61 

CHAPTER S; 

BVALUATION OP 'DIFFERENT ALGORlTHI!^S 

AND CONCLUSIONS 

63 

S.1 

Independent Equation Solution Method 

63 

8.2 

'LO' Algorithm 

63 

8.3 

’PV» Method 

65 

8.4 

PDS and Lumped ¥LS Algorithm 

66 

REFERMCES 


6^8 

APPHDIX 


71 



CHAPTER 1 


INTRODUGTIOH 


State estimation in power systems is a data processing 
algorithm which estimates the true state of the system from 
the following sources of information; 

1 . Measurements of systems variables. 

2. System model. 

3. Apriori knowledge of certain variables known as 
pseudomeasurements. 

The discrepancies in the estimated and true state of 
system are because of meter in-accuracies, communication errors 
errors in pseudomeasurements etc. Such an estimate is necessar 
for deciding various control policies to improve the performan 
of power system in both steady and dynamic states. 

1.1 MOTIVATION 

Many real time control schemes, in fact all of them, 
for static or quasistatic and dynamic performance of power 
systems are benifited by the state estimation. Some of them 
particularly dynamic performance schemes would be impossible 
without state estimation. Some of the real time control schen 
applied in steady or quasisteady state are: 

1 . Steady state security analysis and security 
monitoring. 
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2. Contingency evaluation. 

3. Economic load dispatch. 

4. load frequency control. 

The first tvjo functions cannot be performed in 
absence of state estimation and the latter two functions work 
more reliably with state estimation for the quasisteady state 
being available. The information which is not directly 
contained in the state vector, could be calculated from the 
state vector upon interrogation. 

The state estimation also enables displays ^/rhich 
provides the dispatcher with the following information: 

1. Estimates of lineflows, bus inoection, voltage and 
bus angles etc. throughout the system. 

2. Alarms to indicate system anamolies like information 
system malfunctioning etc. and their geographical loca- 
tion. 

3. Security indices. 

4. Eault location indication. 

State estimate being more complete and accurate, its use for 
data logging is obvious. 

Because of complexities, dynamic state estimation 
generally is confined to estimation of the state of individual 
generators under transient conditions. The need for estimation 
arises because of the fact that the variables of the synchronous 
machine state space models (for example, flux linkages, direct 
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and quadrature axis voltages, damper lar currents etc.)» 
are not directly measurable. The state vector is built by 
measurable output variables of the synchronous machine by 
continuous or discrete time, linear or nonlinear observers. 

The state vector thus built is used to generate an optimum 
governer exciter control to improve the transient response 
of the generator. 

1.2 SUIMARY OP YAxlIOUS CHAPTERS , 

Various methods for state estimation proposed in 
literature amd some modifications in them have been discussed 
and implemented for a standard 14 bus IEEE system in Chapters 2, 
3 and 4. All these methods work on the principle of weighted 
least squares estimation. 

The basi clumped weighted least squares method and 
its fast decoupled version is given in Chapter 2. In addition, 
an estimator based on real power flow and voltage magnitude 
measurement only has been given. Such an esti-iator is 
particularly suited for Indian conditions. Chapter 3 
discusses a very important state estimation algorithm 
which takes only line flow measurements into account. The 
properties of this algorithm which make it suitable for 
online applications are discussed. An approximate method 
for state estimation by solving independent equations is 
proposed in Chapter 4. In spite of being an approximate 
method, it has some plus points over other methods discussed 
in Chapters 2 and 3. 
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Detection and identification of inadvertent bad data 
coming in the algorithm due to conmuni cation failures et», 

IS the subject matter of Chapter 5. The capabilities of the 
detection and identification algorithm are demonstrated on 
the 14 bus IEEE system, Ihe PY method of state estimation 
given in Chapter 2 was used for testing the al-gou-ithm. 

A discrete time extended linear observer to estimate 
the transient state of a synchronous machine given in [ 25] 

IS implemented in Chapter 6. Such observers are necessary 
in cases where a linear observer does not have adequate 
capability to estimate the transient state. 

Chapter 7 briefly discusses the areas related to 
state estimation where a need for further work remains. 
Einally, a comparative study of the static state estimation 
algorithms given in Chapters 2,3 and 4 is made in Chapter 8 
of evaluation and conclusion. 
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CHAPTER 2 

LUMPED WLS AND PAST DECOUPLED STATE ESTIMATOR 

The theory of statistical least squ'TPs estimation 
(lineir estimation) is applied to power systems to produce 
best possible estimate of the true state from the information 
available which is contaminated by metering -^nd communication 
errors. This estimate provides data base for many central 
control and d3 spatch functions as already discussed in 
Introduction. The equations encountered in po^^rer system 
state estimation are nonlinear. They .ire linearized and the 
linear estimation carried out iteratively till convergence 
IS obtained. The computational time and storage requirements 
which are important considerations in any real time control 
scheme , are minimized by ap^-lying decoupling concept to 
st me estimation. 

It is also shown that state estimation can success- 
fully be carried out with real power flow and voltage 
magnitude measurements only. These measuronents are readily 
available hence it is convenient to carry out state estimation 
with them. 

2.1 STATEMEUT OP THE STATE ESTIMATION PROBLM: 

Complete information about the state of a power system 
is contained in the voltage magnitude and bus angles at all 
its buses -with respect to slack bu s. Hence the state vector 
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IS defined as the voltage magnitude at all buses and bus 
angles "with respect to slack bus, Por a ’n’ bus system 
the order of the state vector will thus be ’ 2n-1 ' . The 
state vector in mathematical terms is defined as: 

Vg, . . . , 9 2» ® 3» ®n^ 

where i = 1 , . .n are the voltage rnagnitudes and 
0 i = 2,..n are the bus angles with respect to slack bus. 


In load flow studies also the solution vector is 
same as the state vector used in state estiaation. However, 
numbers of equations is equal to the number of unknown 
quantities and the known quantities In forming the equations 
are assumed to be correctly known. However, in case of state 
estimation, many redundant mea':urements are made so as to be 
able to solve for the state vector oven in the event of 
unavailability of a few important measurements and to facili- 
tate filtering of measurement and communication errors present 
in the available information. The redundancy is defined as 


Redundancy = 


Humber of measured quantities 
Order of the state vector. 


The possible sets of measurements that can be made in 
a system of 'n* buses -ind *1’ lines is given in Table 2.1, 

It IS seen from the table that the maximum figure of 
redundancy is = (6 1 + 4 n)/(2n - 1). It is desired to have 
a low redundancy from economic point of view. However, the 



7 


actual figure of redundancy depends upon technical consi- 
derotions. It has heen found for po’ver systems that a 
figure hctween 1.5 to 2.5 is satisfactory dependent upon 
the system configuration and type of measurements made. 


Table 2.1 


Symbol 

Physical quantity 

Max. number of 
measurements 

^ik'^ik 

Active and reactive power 

flov/s at both ends of the 

lines- 

4 1 


Magnitude of current in 

both ends of the lines. 

2 1 


Magnitude of the bus 
voltage at all buses. 

n 


Reactive and active 

bus injections at all buses. 

2 n 

h 

Magnitude of nodal current. 

n 


2.2 THE LUMPED ESTIMATOR: 

The measurements Z which could be a set or subset 
of quantities listed in Table 2.1 are nonline'r functions of 
the state vector, X, of the form; 

Z = h(X) + noise (2,1) 

The noise is assumed to be zero-mean Gaussian random variable , 
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These equations are linearized at and linear 
estimation theory applied directly to solve them. The 
linearization is done by Taylor series expansion of the 
function h(X) about and truncation upto first order term 
as follows: 

Z = h(X° +AX) “h(X°) + h’(X°)AX (2.2) 

Z - h(X°) = h' (X°) A X 

A Z = H(X°) A X (2.3) 

where AZ = z - h(X°) and H(X°) = h’(X°). 

As seen from equation (A. 6) of J^ppriidn-x , the solution to 
the linear equations (2.3) is given by 

A X = [ H^(X°) ¥ H(X°) ]"^ [ H®(X°) ¥] AZ (2.4) 

The equation (2.4) can be written as 

A X = (2.5) 

where G= H^(X°)¥ H(X°) and b = H^(X°)¥AZ. 

^ IS called information matrix or gain matrix. It is 
symmetric. It has the properties of Gramm matrix when ¥ 

IS taken as identity. 

The equation (2,5) is solved iteratively for AX, 
each time updating the volue of X°, till thel|AX Ik, 
less than a specified limit. 
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2.3 MODEL EOR POWER SYSTEM NEWORK: 

Por simulation tho following types of measurements 
were considered: 

(1) Bus injections, 

(2) line flows, 

(3) Bus voltage magnitudes. 


The models h(X) for power system network which correlates 
the state vector and the above type of measurements are given 
below, 

(a) Bus injections: The network is represented by bus 
admittance model. The bus injections are related to state 
vector by the following expressions. 




n 


Qr, =l ■^r^I I - B^„ COS 0^^) \7^\ 

P P qii P9 P9 P9 P9 9 

p = 1 ,2, , . . ,n 


( 

( 2 . 6 ) 


where P - real power at bus p, 

Jr 

Qp - reactive power at bus p, 

Vp,Y^ - voltage magnitudes at buses p and q, 

Qpq - bus angle 0p - bus angle 0 
^Pl ^ ^Pl” imaginary parts of the pq'th entry in 

'^Bus'*' 

(b) Line flov/s; A line or a off -nominal ratio transformer 
between two nodes, i and k is represented as a four terminal 
network by driving point and transfer admittance (Pig. 2,1). 
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The figure is valid only for line. However it could be 
modified for a transformer. The node voltages ?. and 
are related to currents by 


T 

: 

«i 

hi 

hk 


- 

h 



s 


t 

\ 


The individual terms of the Y matrix are given by 


( 2 . 7 ) 



1 





"ik - \i 


1 



kk 


+ 3 ^ 


( 2 . 8 ) 


The relation between line flows and the voltages at the two 
ends of a line is 


- -* 


= ai... = V- n Y.-+ V. v,. Y 


-T Y. Y 
1 1 11 


1 k ^ik 


( 2 . 9 ) 


The opposite end line flow 


- -* 


- -* 


_ _* _* 


^ki * 3«ki = \Iki = \ h \ Y^^ 


( 2 . 10 ) 


If is expressed as (^ii+0\j_) and Y^^ as (<Ii]j;+OB^jj.) the 
expression for line flows with voltages as polar quantities 
and \Z.®]£ becomes ^ 


= I l" ®11 - I hi l^k I <®ik e ik hk e ik> 

hk = - 1 h f hi " I hi ihi (hk "«\k - hk hk) 

( 2 . 11 ) 


where 9^^ = ( 
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Similar expressions for and 0^^ can le written by 
interchanging i and k in the above expressions. 


(c) Voltage magnitudes: This measurement directly gives 
the bus voltage magnitude. Hence the h(X) is a null matrix 
augmented by identity matrix, 


2.4 ENTRIES OE THE JACOBIAN: 


The measurenent vector can be arranged as, 



^active 


Active bus infections 
Active line flows. 

z = 

^reactive 


Reactive bus in ections. 
Reactive line flows. 


^voltage j 


Voltage magnitudes 


J 1 

1 


The Jacobian for each of these measurements is 
discussed below, 

(a) Bus infection: The Jacobian for these measurement 
IS nothing but the Jacobian used in Ncvjton-Raphson' s method 
to solve load flow problem in polar form with the only 
difference that there is no voltage controlled bus in this 
case. The state vector is taken as [ A9 , |AV| /| 7 |]^ as 
it results into symmetry in certain terms of the Jacobian. 
To summarize the expressions: 


H 

N 

M 

i 

E 


wherg 





15 


for V ^ q. 

“ ""pa = I '"p I P4 ■ V“" ®pq' 

•'pa = -«pq= I ''p H ® PI " ® pa ’ 


for p = q. 


H 

= -Q 

- 3 


PP 

P 

PP 

' P 

li 

s Q 

- 1 

PP 

P 

PP ' 

TT 

= P^ 

+ 1 


PP 

P 

pp 

P 

M 

= P^ 

- ^ n. 1 

V f 

PP 

P 

PP ' 

P ' 


where Pp and Qp are given by (2.6). 


( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


(b) Line flows: The entries of Jacobian corresponding 
to this measurement are found by taking partial derivatives 
of the real and reactive line flows of equation (2.11) with 
respect to bus angles and volta^ e magnitudes. 1tfe get eight 
terms viz. 


3 P 


ik 


3 0 1 


36 


3 Q 


ik 


3 0 


3Q 


9«lk II 

3 I V^l ’ 3 ® 


ik 




ik 


TIvTT ' 


These derivatives arc then to be multiplied by ] 7 | 
wherever necessaiy because the state vector is assumed to 
be (A0 » a| 1/ 1^ I )^* performing the above operations, 


the follo\ving expressions arc obtained: 
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= sin 6^^ + oos 


(2.18) 


2|"il Sii-l^ill\|CSik“^\k + B.^=-eik) (2-19) 


'^I'V ^9lk '•ik + 2^k 


(2.20) 


'"i' "i' ' ■"‘^' '°xk^“\k-\k“^ ik^ (2.21 ) 


’^Ik 


.S^lk. 


iTj|vJ(a^k ^"“«lk-Blk (2.22) 


(g-jfT^I^I = I \l I \l (9ik 9 ik + Sik 9 ik) 

Ic 


(2.23) 


9 Q-, 

^ik 


" I ^il i ® ik “ \k ® ik^ (2.24) 


(g-j^^l V = I ’^ll ('9lk =1-n9 oos 6.^). 


(c) Voltage magnitufTes: The partial (ilerivatives of the 
function relating voltage magnitude mea 'urements and the 
voltage magnitudes, -with respect to V vill he unity and 
v/ith respect to 9 is taken as zero. Thus the Jacobian for 
this mea urement becomes: 

(2.25) 


Null matrix 


Diagonal matrix v/ith 
diagonal entries 
j equal t o 

I voltage magnitudes 
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To summarize the overall Jacobian for the three kinds 
of measurements discussed above: 


load flow 

Load flow 

H 


Line flow 

Line flow 

1 D P , 

ik . ik 

3P , 3 P , 

( 1 |y I '( Iv 1 


3 Y. ^ 9 V J ^ 

1 k 

Load flow 

Load flow 

M 

I. 

Line flow 

1 

1 Line flow 

" hk . 8 «ik 

(\Ak)l7 1 

9 01 ’ 3 Qic 

iU.\ ^ 

X K 

Zero 

Diagonal matrix with 

j 

entries Y^ in the 

' 1 

diagonal 


(2.26) 

2.5 IMPIMEFTATION: 

The lumped ^VLS method was implemented on the IEEE 

14 bus test system. The results obtained from load flow 
analysis were used as input data. The weighting matrix was 
chosen to be identity. The flow-chart for implementation 

15 given on next page. Results obtained are listed in 
Table 2.2. 
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FLOW OHiRT FOR lUMPEl? ¥13 METHOD: 


Take flat start for the state vector "f 

6=0, V = 1 .0 

j 

— . — 

Specify convergence limits e-] and £2 
for angles and voltages. 

—I 

Calculate h(X°) 

f ^ 

Build Jacohian H 


Evaluate (H^) and H^(Z-h(X°)). 
let G = and h = H^(Z-h(3[P)) 


T 


Solve 

G Ax = h for Ax hy triangular 

factorization 


I 

1 







Update X = X+ Ax 




DE1TA1 = Halloo 

DE1TA2 = If V|I 
' *00 
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Table 2,2 

RESULTS OE THE lUMPEE WIS METHOD 


BUS 

HO. 

LOAD EL0¥ 

SOLUTIOH 

VLS 

SOLUTIOH 

Yoltage 

tingle 

Yoltage 

Angle 

1* 

1 .06 

0.0 

1.06 

0.0 

2 

1 .045 

-4.93 

. 1.03155 

-4.65657 

3 

1 .010 

-■»2.72 

0.99398 

-12.34901 

4 

1 .019 

-10.33 

1.00330 

-9.96780 

5 

1 .020 

-8,78 

1.00427 

-8.37051 

6 

1 .070 

-14.22 

1.05344 

-13.91745 

7 

1.062 

-13.37 

1.05811 

-13.14338 

8 

1 .090 

-13.36 

1.08300 

-13.12710 

9 

1.056 

-14.94 

1.05811 

-14.75851 

10 

1 .051 

-15,10 

1.05586 

-14.99860 

11 

1.057 

-14.79 

1 .05310 

-14.75348 

12 

1.055 

-15.07 

1 .05524 

-15.13332 

13 

1.050 

-15.16 

1 .04792 

-15.16305 

14 

1.036 

-16.04 

1.04362 

-16.09178 


Redundancy = 3,0 

Ho. of iterations = 4 

Total time taken =133 sec. 

= 0.0005 

£2 = 0 . 0001 = 

* indicates slack bus. 
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2.6 THE EAST DECOUPLED ESTIMATOR (EDS): 

The concept of P-0, Q-V decoupling used in load flow 
studies IS extended to state estimation. The hasic work: in 
this area is done in ref . [2 7]. The information matrix Or of 
eq.(2.5) is approximated ghy two suhmatrices each correspon- 
ding to the respective decoupled loops. This reduces the 
storage and also the total time taken to converge, though 
the total number of iterations are more than that required by 

the lumped WLS estimator. The reduction of computational 
* 

time IS due to the fact that the submatrices derived from the 
gain matrix G are treated as constants. They need not be 
triangularized every iteration. 


The measurement vector Z arranged in active and 
reactive groups is expressed in termsof state vector as 


The Jacobian H(X) of the function h(X) is given by 



= h(e , V) 


(2.27) 


active 


reactive 


§h ^ 3 h 3 h 

dX "SF" Ty 


3f/ 90 

3 f/ 37 

9g/ 9 6 

3 g/9 V 

I 


^11 

■ ^12 

H21 ^ 

CM 

CM 

W 

I 


(2.28) 


The information matrix G of eq.(2,5) is 


R-'’ ^ r~''-2s 

90 a 90 90 ^ 9© 

m!r-i » + l^R-i 

0^ a 3 ^ 90 r 9T 

R-’' + i£ R-1 1_S 

37 a 30 ^ 37 r 30 

3f^ R-1 3f ^ 3g^p-1 9£ 

37 a 97 97 ^r 97 


77 (2.29) t 


G(e,V) = 
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vhere and R^^ are diagonal matrices and 


= R = 


R 

a 

0 

0 



W is the weighting matrix. 


She above information matrix G( 0 , "V) is approximated by 


G( ®, V) = 


3f® j^-1 3f 

30 a 3Q 

0 

0 

37 r 37 


(2.30a) 


by assuming that the matrices and H 22 are large compared 
to 2 ^ 21 ' eqn.(2.28) for information matrix only. 


A further simplification in the information matrix 

computation is made by calculating it once at the previously 
value of j 

estimated /_ X or flat starting value of X, i.e. V=(1,1...1) 

T 

and 0 = ( 0 , 0 , ... 0 ) , The information matrix then becomes 
constant independent of 0 and V. 


j ^ 

0 


0 I 

B 



(2.30b) 


Thus the iteration scheme becomes 


r 


1 


0 


(AX) = R“\z-h( 0,7)) 


(2.31) 


[0 B ] 

or A( A 0 ) = ,7), e,V)] R"^(Z-h(e,V)) (2.32) 

B(A7) H2^0,7)] R“\z-h(0,7)) (2.33) 

The equations (2,32) and (2,33) can be solved iteratively to 


get the optimal solution of state vector X, The two matrices A 
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and B are symmetric and sparse in nature which need he 
factorized only once, 

2.7 IMPIMMmTION; 

The PDS algorithm was implemented on the IEEE 14 bus 
test system. Results are listed in Table 2.3 along with the 
true solution obtained by load flow. The flow chart for 

I 

implementation is given is given below. 

FLOW CHART FOR EDS ALGORITHM IMPLEMENTATION: 
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0?a]3le 2.3 

RESULTS OP EDS METHOD 


Bus 

No. 

load flow solution 

EDS solution 

Voltage 

Angle 

Voltage 

Angle 

1*" 

1 .06 

0.0 

1.06 

0.0 

2 

1 .045 

-4.98 

1,0316 

-4.6561 

3 

1 .010 

-12,72 

0.9942 

-12.3448 

4 

1 .019 

-10.33 

1 .0036 

-9.9639 

5 

1 .020 

-8.73 

1.00456 

-8.3675 

6 

1 .070 

-14.22 

1.0558 

-13.9032 

7 

1 .062 

-13.37 

1 .0536 

-13.1329 

8 

1.090 

-13.36 

1.0335 

-13.1161 

9 

1.056 

-14.94 

1.0586 

-14.7456 

10 

1.051 

-15.10 

1.0564 

-14.9851 

11 

1.057 

-14.79 

1.05362 

-14.7393 

12 

1.055 

-15.07 

1.0556 

-15.1654 

13 

1 .050 

-15.16 

1.08484 ' 

-15.1468 

14 

1 .036 

-16.04 

1 .04419 

-1647693 

Redundancy 

= 3.0 



Uo. 

, of iterations 

= 12 



Total tide taken 

=a 99.66 sec. 



El 

= 0.002 




^2 

= 0.001 




* indicates slack bus. 
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2.8 ACCEIiERATIOU OP OON "ERGENOE: 

The above approximations are advantageous in terms 
of storage and also the total time required for convergence. 
The only short coming is that it takes relatively more 
number of iterations to converge. Attempts have 
been made co minimize the number of iterations for such 
a fast decoupled method [2?]. The scheme suggee ted in [ 27] 
is as follows: 

The iterative scheme is written as 


h+1 

•r-i 

CD 

II 

* ^ 

(2.34) 


= 7. 

1 

+ T. 

1 

(2.35) 

where AS^^ 

= 

^ ( 9, , T ) 

39 ^ ^ 

(2.36) 

B 


( 0,, V ) 

37 11 

(2.37) 


The convergence is accelerated 

size u. and 3. as^follows: 

1 1 






5 = e. + a s 

1+1 1 11 


(2.33) 


Y .= V. + 
1+1 1 


3. T . 
1 1 


(2.39) 


The choice of and 3^ should be theoretically made 
su ch tha t 


J( 0^ +a ^ S^.V^) = min {j(0 +« , 7^) } 

and 


J(7^ + ^ + ^'^i 

3>.0 


Vi)> 


(2.40) 


(2.41) 
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and can be found out by a linear search or curve 
fitting. A simple scheme to choose and 3^ can be 
setting 


1+1 


(2.42) 

and 

a . . 

1+1 

= a./2 If V^) > J( 0^, V^) 

(2.43) 


Similar scheme can be used for 3^. This simple scheme is 
reported to have reduced the number of iteration by 20^ 
in[27] • -ftficexFtiy oohome -feae "ispeesi s uggest ed 4s 


2 . 9 STATE ESTIMATION BASED ON READ PO'^i^ER ELOW IN LINES 
AND VOLTAGE MAGNITUDES ONLY (PV) : 

State estimation can successfully be carried out 
with real power flows in lines and voltage magnitude measure- 
ments only. These kinds of measurements are readily available 
and state estimator based on these measurements will parti- 
cularly be suited for Indian conditions. 


The derivation of algorithm directly follows from 
eqn.(A.6 ) of Appeodrs . The Jacobiansfor these kinds of 
measurements are given in Section 2.4b and 2.4c. The over- 
all Jacobian will have the following form: 


H = 


Line 

flow 

3P , 

3P , 

ik 

ik 

36 , 

3'Ot 

1 

k 

Null 


Line flow 

(_^) V. , (-^,) vj 
3| V I 1 thj ^ 

X K! 


Diagonal matrix with 
entries 1 V 1 on the 
diagonal 
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Iterative scheme given in eqn.(2.4) is then used to 
solve the state vector. 

2.10 IMPLEMOTATIOF: 

The PV state estimation was carried out for the 
IEEE 14 bus system. Results are tabulated in Table 2.4. 
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Table 2,4 


RiSUlTS OP PV MITHOD 


Bus Ro, 

load flow solution 

PYO solution 

Yoltage 

Angle 

Yoltage 

Angle 

1* 

1,06 

0.0 

1 .06 

0.0 

2 

1.045 

-4.98 

1.0437 

-4.965 

5 

1 ,010 

-12.72 

1 .0105 

-12.727 

4 

1 .019 

-10.33 

1 .0217 

-10.398 

5 

1 .020 

-8.78 

1 .0180 

-8.762 

6 

1.070 

-14.22 

1 .0782 

-14.932 

7 

1.062 

-13.37 

1.0582 

-13.621 

8 

1 .090 

-13,36 

1 .090 

-15.621 

9 

1.056 

-14.94 

1 .0506 

-15.227 

10 

1 .051 

-15.10 

1.0503 

-15.488 

11 

1.057 

-14.79 

1.0565 

-15.229 

12 

1.055 

-15.07 

1 .055 

-15.563 

13 

1.050 

-15.16 

1 .0498 

-15.610 

14 

1 .036 

-16.04 

1.0359 

-16.495 


Redundancy 

Ro. of iterations 

Total time taken 

^2 


= 1,27 

=r 4 

= 86, 6S sec, 

= 0.0005 

= 0,0001 


* indicates slack bus 
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Table 2.4 


RiSULTS OF PV MITHOD 


Bus No. 

load flow solution 

PVO solution 

Voltage 

Angle 

Voltage 

Angle 

1* 

1.06 

0.0 

1.06 

0.0 

2 

1 .045 

-4.93 

1.0437 

-4.965 

3 

1 .010 

-12.72 

1 .0105 

-12.727 

4 

1.019 

-10.33 

1 .0217 

-10.398 

5 

1 .020 

-8.78 

1 .0180 

-8.762 

6 

1.070 

-14.22 

1 .0782 

-14.952 

7 

1.062 

-13.37 

1 .0582 

-13.621 

8 

1 .090 

-13.36 

1.090 

-13.621 

9 

1 .056 

-14.94 

1 .0506 

-15.227 

10 

1 .051 

-15.10 

1.0503 

-15.483 

11 

1.057 

-14.79 

1.0565 

-15.229 

12 

1.055 

-15.07 

1 .055 

-15.563 

13 

1 .050 

-15.16 

1 .0498 

-15.610 

14 

1 .036 

-16.04 

1.0359 

-16.495 


Redundancy 

No. of iterations 

Total time taken 

^2 


1.27 

4 

86. 6S sec. 
0.0005 
0.0001 


* indicates slack bus. 
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OHAPIER 3 

SHE LINE OKLY ALGORITHM (’LO' AL&ORITHM) 


A very efficient state esti.iiation algorithm has hern 
developed in [34j‘which takes into account only line flow 
measurements. It is superior to lumped WLS and EDS algorithm 
in the following respects: 

(1) Easter convergence. 

(2) Simplicity of the algorithm and ease with which 
changes in network configuration or metering 
location can he taken into account. 

(3) Lesser storage requirements. 

However this algorithm has lesser flexibility for 
taking into account differing meter accuracies compared to 
lumped ¥1S and EDS algorithms, ii- further sacrifice in the 
flexibility to incorporate differing meter accuracies pays 
in terms of reduced storage when line flows in both ends of 
all measure'^ lines are available. 

3.1 THE AEP LINE OHIY ALGORITHM 

A line between buses i and k is represented by a 
four terminal equivalent model as given in Section 2.5b 
of Chapter 2, Equations (2.7), (2.9) and (2,10) are 
repeated here for easy reference. 



27 





^ii 



V. 

1 






* 




1 



k. j 

-* 

^ik"-^^ik = 

^ I 
i -^ik 

= Y. 

1 

7 . Y. . 
1 11 


Y Y Y 
1 k ^ik 


- -* 




= \ ^ki = \ ^ \ \ ^kk 


(2.7) 

(2.9) 

( 2 . 10 ) 


The terms have already been defined in Sec. 2. 5 of Chapter 2. 
The measured quantities are real and reactive line flov/s 
from bus 1 to k or k to i or both. These measured quantities 
are translated into voltage difference between the two nodes 
1 and k. line flow translated into 3-nd line 

flow the following relations: 

- - - I’.m-dQ 


TT Y — Tj ik li ik ik 

V -V, - V. — 3 -^ 

X. X Y V Y 

ik i ^ik 


(3.1) 


and 


- Y, +Y, , P, .-oQ, . 
"vr ir ^T ki kk ki ki 
Vi.“V_, = Vi, ::: ry ~r » — 


k 'i 


(3.2) 


'ki 


\\i 


These translated quantities are treated as measured 
quantities in further derivation. Call the complex vector of 
the translated measurements as The function h(X) which 

correlates and X is 


_-h(X) = B^X - B 27 . 


R 


(3.3) 


where B^ is the bus incidence matrix and B 2 is the column of 
element node incidence matrix which is deleted to obtain bus 
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incidence matrix B-j and is the complex slack bus voltage 
assumed to be correctly kncwn. 

To apply equation (A,6 ) of ilpptsnd'ix to estimate the 
state vector from the above translated measurements the 
Jacobian H(X) is calculated by differentiating equation (3.3) 
with respect to X ; 

H = (3.4) 

The Jacobian is constant. Thus the state estimate is obtained 
by iteratively applying equation (■^•6 ) of Appendix as follows 

AX = (B^ ¥ B^)"^ (B^ W) (U^ - h(X)) 

= (B^ W B^)“'' (B® ¥) (U^ - B^X + \) (3.5) 

where ¥ is the weighting matrix for measurements. It should 
be noted that equal weights are assigned to P &*Q measurements 
corresponding to individual line. 

Equation (3.5) is simplified by expanding the expression 
as 

A X = (B^ ¥ B^ (B^ ¥) (U^ + B 2 \) - X 
or 

as X IS X +A X. 

new 

The equation (3.6) is solved iterativ ely replacing T 
by until convergence is obtained. 
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3.2 M03DIFIC-.i:i01I OP AEP iiLGORi™ TO TiJCE IFTO ACCOUNT 
EPPICIENTLY BOTH EtTU LINE EDOWS 

When information about both end line flows in all 
the measured lines are available the algorithm of Section 5.1 
can be modified which will reduce the size of matrix to 
half of what will normally be required by equation (3.6), 

The assuraption made here is that weighting factor for the 
four measurements, viz, P and Q at both ends of a line are 
equal. This approximation is worth considering as it 
results in significant reduction in storage. 

The derivation is as follows: 

Let be the vector of complex element voltages 
translated from measurement using equation (2.10) 

and Uj^ as described earlier be calculated from measurements 
Pij^+jQi^ using equation (2.9). We define two functions h^(x) 
and h 2 (Y) similar to equation (3.3) relating state vector 


X with and U^ respectively; 

h^ (X) = X - Bg \ (3.7) 

and hgCX) =-B^ X - ^2 \ 

The Jaconian for these functions are 

= B-j (3 .9) 

H 2 = (3.10) 


Equation ( A. 6) of Appendixi is modified to take into account 
n different types of measurements as 



where - Jacohian corresponding to ith type of measurement, 

- Weighting matrix corresponding to ith type of 
measurement 2^, 

h^(X)- Model for ith type of measurement Z^. 

Considering and to he two different types of 
measurements with equal weighting matrices we get, 

X = )“"* [ B^(U^-B^X+B2fj^)-B^¥(U2+B^I+B27j^) ] 

= (2B^^re^)"^ [ B^W(U^-U^)-2B^B-^X + b3[w(2B2\)1 

= (B^¥B^ )“^ { ¥ B2 \ ) - ^ 

\ew = (=1 '' ’■'' “l W(^) H- Bj Tjj 1 (3 -11) 

Equation (5.11) is solved iteratively till convergence is 
obtained. 

In order to have a complete state vector, it is 
necessary to have metf~rs placed at least on lines that form 
tree of the v;hole netv'ork. When algorithm of equation (5.11) 
is being used, we have to have both end line flows measured. 

In case one end line flov; is missing, the missing measurement 
IS put e lual to the other end line flow available. Thus the 
minimum figure of redundancy for the algorithm of equation(5.H) 
IS 4n/2(n-1) (approximately 2.0) for an n bus 1 line system. 
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It IS possible to apply this method to parts of the 
network if each such part has one reference bus. 


3.3 liiPIiEMMIATION 


The 'LO* algorithm discussed above are very simple 

I! 

to be programmed. The gain matrix can be built just 

by inspecting the net ring configuration as is built. 

It IS very sparse. Its structure is same as that of of the 

metered part .of the network. Off-diagonal entries are either 


1.0 or -1.0 when is taken as identity. The following 
algorithm builds the gain matrix, when W is identity, following 
is the description of variables used. 

(a) Bor an n+1 bus and 1 line system starting end index of 

lines are stored in vector NX(I) and terninating end index in 
NY(I) for I = 1 to 1. (n+1)th bus is slack bus. 

(b) The product (B.|B.j) is stored in B. 

T ( 

Algorithm for finding (B>|B.j): 


Step 1; Initialize B(i, 3) = 0, i = 1,...n, o = 1,...n, 
set JOC = 0, i = 1 , 


Step 2: Set KX = KX+1 , KY=0. 

Step 3: If NX(i) = KX go to step 4 else if NY(i)= KX 

go to step 6 else go to step 8. 

Step 4: KY = KY+1 , n1 = ]^Y(i). 

Step 5: If n1 n+1 then B(KX,n1 ) =-1.0 > 

Step 6: KY = KY+1, n1 = WX(i). 

Step 7: If n1 n+1 then B(KX,n1 ) = -1,0 » 

go to step S. 


go to step S, 



32 


Step 8: i = i+1 . 

Step 9: If i.LB.n go to step 3. 

Step 10: B(KX, KX) = float(KY), 

Step 11: If KX.LE.n go to step 2. 

Step 1 2: Stop. 

T 

Ibe matrix also need not be stored explicitly 
as sufficient information is contained in arrays FX and NY 
to calculate the product of B^ with an array. 

Elow chart for implementation of the ’10’ method is 
given on next page. 
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START^ 


Read the system information and 
measurement vector Z consisting 
of line flows 

Assume a state vector X°,set jiL^= 0 
Rormlate matrix, store in ^ 


Triangularly factorize A 


1 


M5 


Formulate bus incidence matrices 
and vector B 2 


m 

Calculate (B^W) and store in B^ 
itself 


i 


2[K=KK+1 


Calculate and from eq.(3.1) 
and (3.2) 


Find the sum 

store it as a complex vector in 


Find the product B^ and store 
in a complex vector U 


Solve equation of the form aX=TJ. A is 
stored in triangularly factorized form. 
Store solution in itself. 


Find II cabs(X-U) II go = DELTA 
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3.4 COMPUTER RESULTS 

The algorithm (5.3) was implemented on the IEEE 14 bus 
test system. The results obtained and details of the algorithm 
are listed in Table 3.1. 

Table 3.1 


Bus 

Wo. 

load flow solution 

'LO' 

algorithm 

Voltage 

Angle 

Voltage 

Angle 

1 * 

1.06 

0.0 

1.06 

0.0 

2 

1 .045 

-4. 93 

1.043 

-4.93 

3 

1.010 

-12.72 

1 .007 

-12.757 

4 

1 .019 

-10.33 

1.015 

-10.356 

5 

1.020 

- 8.73 

1.022 

- 3.770 

6 

1 .070 

-14.22 

1.002 

-14.922 

7 

1 .062 

-13.37 

1 .029 

-13.6*28 

S 

1 .090 

-13.36 

1 .053 

-13.628 

9 

1.056 

-14.94 

1 .014 

-15.402 

10 

1 .051 

-15.10 

1.009 

-15.570 

11 

1.057 

-14.79 

1.015 

-15.244 

12 

1.055 

-15.07 

0.993 

i -15 .303 

13 

1.05 

-15.16 

0.995 

-15.309 

14 

1 .036 

-16.04 

0.936 

-16.689 


Redundancy = 2.61 

Total Uo.of iteration = 5 

Time taken = 19.7 sec. 

e = 0.00009 


* indicates slack bus. 
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CHAPTER 4 

STATE ESTIMATION BY IHDEPEHDEHT EQUATION SOLUTION 

A method is proposed for state estimation which 
makes hus-hy-hus sequential estimate of the state vector 
thereby obviating large matrix storages. Also this method 
IS faster than the methods discussed in Chapters 2 and 3. 

The bus-by-bus sequential estimate is obtained by splitting 
the lumped equations into independent equations. Accuracy 
of such an independent equation solution will depend upon 
network topology and instrument location. Such a technique 
when applied to a radial system may give unacceptable results 
as the error in previous bus voltage estimate will go on 
propagating. However, for systems with "favourable geometry" 
"good local redundancy" at all the buses and"properly chosen 
instrument location" reasonably good estimates can be 
obtained. 

4.1 SPLITTIHG lUMPEE EQUATIONS INTO INEEPEMEENT EQUATIONS 
The line flows at both ends and voltage magnitudes 
at all the buses are considered to be measured quantities. 
Like in the 'LO’ method, here also the slack bus voltage 
IS assumed to be accurately known. The measurement vector Z 
IS arranged as follows: 

(1) The line flows at both the ends of all lines 
connected to slack bus and the voltage magnitudes of the 
buses on which these lines terminate come first. 
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(2) Next similar sets of measurements come, for lines 
connected to the buses on which the previous lines terminate. 
This is carried out till all the measured lines are considered. 

Consider as an example the part of network shown in 
Pig. 4.1. The measurement vector Z consists of line flows 
S-j,S 2 ... etc. and voltage magnitudes 7^,72*.. etc. This vector 
is rearranged in the following way: 

= (S^, S2, 7 ^ : 83,3^, 73 : S3, Sg, 7^ : ...) 

( 3^ . ^2 ; ^3 . » . . ) 

The function 3 = h(X) takes the following form. 

z, = h(X,, V^) 

Zj = hCXj, v^) 

2^ = h(X2> ) 

where 7^ is the known complex slack bus voltage and X.^,X 2 *-* 
are the complex bus voltages of buses 1,2 ... . 

Each of these measurements are independently used to 
estimate one bus voltage which is treated as constant in 
subsequent bus voltage estimation. In mathematical terms, 
this amounts to mean that we obtain X..j,X 2 ,.... sequentially, 
such that they minimize the local cost functions J..^,J 2 >J^ etc. 


(4.1) 


(4.2) 


defined as 



= (Z^-h(X^,V^))^ ¥(Z^-h(jC^,V^)) 
J 2 = (Z2~h(X3, Y^)f ¥(Z2-1i(X3,V^)) 
J3 = (Z^-h(X2,X^))^ ¥(Z3-h(X2,-S^)) 
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( 4 . 3 ) 


Strictly speaking *^1 and ^2 and and ♦.. etc. 

min min 

does not mean J(X) (the total cost function defined in WjS 
min 

■nethod). However, since the local redundancy for each 
measurement set is high (2.5) results will not deviate 
much from the VIS solution. 

4.2 SOIUTIOH PHOOEDUKE 

The group of measurements Z^ is considered first. 
The two line flows in both ends of the line connected to 
slack bus to which the group Z^ corresponds, give t^vo 
solutions for the voltage of bus at the other end of the 
line. Third solution to the voltage magnitude at that bus 
IS directly obtained by the measurement. Veighted average 
IS then taken of the three solutions to voltage magnitude 
and two solutions to phase angle. Assuming this solution 
to be final one, the solution for the bus voltage of the 
bus connected to this bus is next carried out in a similar 
way. Procedure is repeated till all the bus voltages have 


been obtained. 
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It IS clear that the deviation of such a solution from 
the solution becomes more and more enhanced as the number 
of measured lines in series go on increasing. Thus a favourable 
geometry will be one where we have as many lines directly 
connected to slack bus as possible. An extreme case which 
of course, is not practicable, is one having all the buses 
directly connected to the slack bus (Fig, 4. 2), For such a system 
the soluiiion obtained by the above method and the WLS method 
(LO) will be exactly same. Results in both the cases are 
sensitive to the error in slack bus voltage. 


In Figs. 4, 3 and 4.4, Fig. 4. 3 will give results more 
close to the WLS results than will Fig. 4. 4 give. The 
importamoe of meter location in obtaining better results is 
emphasised. The instrument location of Fig. 4. 5 is a better 
choice than that of Fig. 4. 6 for the same system. 


4.3 SOLUTION FOR VOLTAGE FROM LIWF FLOWS 

The measured line flows from bus ’i' to 'k' and 
'k' to are related with the voltages at buses ’i» and 'k' 










ik 


^ki-^Qki = v^ki " Vk4 


(2.9) 

( 2 . 10 ) 


ki ki 

where I , ••• \k Vk 


meanings. 



a 


IS unknown in iDoth these equations. Eqn.(2.9) 

IS linear when voltages are represented in rectangular form. 
Eqn.(2.10) is nonlinear and has to he solved hy ll-R method. 
In rectangular form the solution to will he 



®k(new) 

(4.5) 

^k(new) ^ ^^k 


where H = 2e^ 

(4.6) 

" = 2fk \k ^ ^ ®kih 

(4.7) 

“ ^®k®kk " ®i®ki ®ki^i 

(4.S) 

!■' = -2fk\k - \rh + \®kr 

(4.9) 


The following terms have been used in expressions for 
H,U,M and I, 


Y. . = 

G.. + 

D®- ,■ 

11 

11 

11 

hk - 

®rk - 

^®ik 

II 

•H 

®kr - 


Rk “ 

G. , + 

kk 

^\k 

h = 



\ = 

®k + 

3^k 


( 4 . 10 ) 


(4.11) 
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As eqn.(2.9) is linear, solution to is obtained 
by solving this equation and this solution is used as 
initial guess for the iterative scheme (4.5) for eqn.(2.10). 

Solution can also be obtained for in polar form. 
Ihe Jacobian for F-R solution in polar form is given in 
Section 2.4b. The advantape of solving in rectangular 
form IS that we get one linear equation >iiioh is easy to 
be solved and since there are no trigonometric functions 
involved any where in the Jacobian, convergence is faster 
than the polar form of solution. Hardware implementation 
of the rectangular version will be economic compared to 
polar version as only basic functions of addition, subtrac- 
tion, multiplication and division are involved in it, 
whereas, additional software implementation of summation of 
series or hardware implementation of PROMs or EPROMs is 
required for SIR and COS functions in the polar form of 
solution. 

# 

4.4 COMPUTER RESULTS 

The IEEE 14 bus test system was used to test this 
method. The results are listed in Table 
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(Tabl e 4 » 1 

EESULTS OP IlTDEPEtrDEfTT EQUATION SOLUTION 


Bus No. 

Load flow 

solution 

Independent equation 
solution 


Voltage 

Angle 

Voltage 

Angle 

1* 

1.06 

0.0 

1 .06 

0.0 

2 

1.045 

-4. 93 

1 .045 

-4.977 

3 

1.01 

- 12.72 

1.00994 

-12.7177 

4 

1.019 

-10.33 

1.02037 

-8.7784 

5 

1.02 

-3.78 

1 *00523 

--7*9923 

6 

1.07 

-14.22 

1*04252 

-20.2528 

7 

1 .062 

-13.37 

1)*C503 

*8.8496 

8 

1.09 

-13.36 

0*9S937 

-15.2685 

9 

1.056 

-14.94 

1.02575 

-13.68259 

10 

1.051 

-15.10 

0.9991 

-19.22794 

11 

1.057 

-14.79 

1.03722 

- 9.442 

12 

1.055 

-15.07 

1.07872 

- 8.85028 

13 

1.05 

-15.16 

1.03045 

- 9.81833 

14 

1.036 

-16.04 

1.01512 

-8 .94556 

Total time 

taken : 2.466 

sec. 



Iterations 

taken : Max. 

3 iteration 

per line. 
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GHAPTER 5 

EEPEOTIOR ARE lEEUTIEIOATIOR OE GROSS MEASURIMMg ERRORS 

State estimation algorithms produce reliable estimate 
when there is no gross measurement error in the input data, 
i^etection and identification of inadvertant had data coming 
in the input data is of great practical importance to enhance 
operator confidence in the state estimates. A sound acceptance 
criterion for results required is' obtained by statistical 
analysis of estimates, 

5.1 OE THE PROBLEM 

In power system model, 

Z = h(X) + 

Vg, IS assumed to be zero mean G-aussian random variable, 

However, a sudden gross error appearing in the input data 
violates the assumptions of Gaussian distribution of v^. 

The model in such a case is written as 

Z = h(X) + V 
where 

V = V + ae ; e = (0, 0, ... 1 ... 0). 

2i 

rp 

Entry 1 in e corresponds to bad data location or locations 
and G IS the size of bad data. The problem is to detect the 
presence of bad data and identify which measurement is in 


gross error 
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Detection : 

Inspection of the cost function J(1t) of the least 
squares estimation, at optimal solution X reveals whether 
the input data is in error or not, as had data is expected 
to cause an unexpectedly large J(X). The question of 
selecting the threshold value of J(X) is answered hy 
statistical properties of J(X). 

Identification : 

Having detected the presence of had data, identifying 
which measurement is in error is the next question to he 
answered. Unfortunately it is not always possible to infer 
correctly about such had data location from the residue 
(Z-h(X)). The largest residue (S-h(X)) need not necessarily 
coincide with the had data location. This phenomenon is 

A 

called smearing. The residual vector (Z*-h(X)) is normalized 
hy variances of its respective components so as to "see 
through" the smeared residues. The variances of individual 
entries in the residual vector are obtained hy stastical 
analysis. 

5 . 2 MTHMATIGAl MCKGHOUHD 

The estimate X is based on the criterion of 
minimizing the cost function: 

J(X) = (Z-h(X))® R“\z-h(X)) 

. = H^R"”' [Z-h(X) ] = H^R“^r =0 (5.1) 

x=x 


'•^hus 


dj 

dX 
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A 

Where H is the Jacobian and r is the residual vector, 
r can be written as 

r = h(X) + v-h(X) (5.2) 

where X is the true state vector. Defining state estimation 
error as 

6X = X - X (5.5) 

It follows from (5.3) and (5.1) 

H®R"\h(X) + V - h(X)) = 0 

Substituting h(X) = h(X) + H6X; 

(H 6X + v) = 0 
or 

6X = -(H^R"^H)"^ H^R"''v (5.4) 

or 

6X = - h'^R'^v with = (H^R“^H)"^ (5.5) 

with V = i.e. zero mean G-aussian random vector; the 
covariance matrix of estimation error 6X is 

E[6x 6^] = (5.6) 

The residual vector r can be written as 

r = h(X) + V - h(X) = W V (5.7) 

where 

¥ = (I - El R”^) 

IS called residual sensitivity matrix. 


(5.3) 
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The residual covariance matrix becomes 

e[ t ^ = {r ^ E H®) = ¥ R (5.9) 

Ihe difference between true and estimated measurements is 
given by 

6 h = h(l) - h(X) = (¥-I)v (5.10) 

¥e now define weighted and normalized residuals as 

^ ^ (5.11) 

where D is a diagonal matrix obtained by 


H = diag 


( 5 . 12 ) 


The inequality r^^ ^ r^ holds. Por v = 


v^ the covariances 


of r-^ and fj^ are 






and 


E Ir^ = Ujj =fD-'' VD-^ 


(5.13) 


(5.14) 


The diagonal elements of % and Ujj. are equal to one. 
Define the normalized residual sensitivity matrix ¥j^ as 


¥j^ = fD""* ¥ 


(5.15) 


If ¥„ , is the element corresponding to gth row and kth 

column^'^¥jj- then 


¥. 


I 


'r.jte 




r,kk ^2 


(5.16) 
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where is the standard deviation of the 3th measurement 
J 

error. Since is positive semidefinite , it follows from 
( 5 * 16 ) that 

^^^,33*1 i f\,3k:l k=l, 2 ,...,m 

where m is the size of the measurement vector. 

This means that the largest nomalized residue r^^ 
coiacides with the had data location. This is not always 
true with 

Considering the function J(X), it is easy to see that 
it has chi-square distribution with k degrees of freedom 
when V = v^, 5 k = size of the measurement vector(m)- size 
of the state vector(n). 

In the presence of had data 


V = v^ + e (X 

the residuals r and normalized residuals take the form 
r = ¥ e a 

and (5.17) 


% = % 


Thus J(X) becomes 


e a 


J('i') = v^ W V + 2a e^R“^ v + e^ R”^ e (5. 18) 
z z z 

The first term is chi-square distributed, second is normally 
distributed and third is a constant. 
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Detection and Identification Theory ; 

Detection and identification is carried out by 
hypothesis testing. Two hypothesis HO and H1 are defined. 

HO : Ho bad data 

H1 : HO IS not true. 

(1) The"J(X) Test": 

Accept HO if JvX) < Y (Y = chi-sq.uare distribution with 

k degrees of freedom and specified 
false alarm probability b defined 
as the probability of making a wrong 
decision of rejecting HO when 
HO IS true). 


Reject HO, otherwise. 

(2) The "rj^ test": 

Accept HO if I rj^ ^ I < y 

Reject HO, otherwise. 

(3) The "r^ test": 

Accept HO if | r^^ k^ ^ 

Reject HO, otherwise. 

The same threshold y being used, r^^ test is more 
sensitive to errors since r^^ always. 

The following recommendations are made for implementa- 
tion of bad data detection and identification tests. 



50 


(1) Use J(X) test to detect the presence of had data. 

(2) If had data is detected, apply or rj^ test to 
identify the had data, r^^ test is more accurate hut 
needs off-line calculation of matrix D. 

(3) -A-pply J(X) test to detect the validity of prevaous state 
estimate at a given instant of time. Use the new 
measurements and the old estiaate X to find J(X). 

If J(X) IS found to he significantly different from 
the previous J(X) a new state estimation is required. 

A plot of such a function J(X) with time may he useful 
for the system behaviour detection. A gradual increase 
in J(X) will mean that the system state is steadily 
changing whereas a jump in J(X) will indicate a 
line outage or a %nir.den' iQa-d ov^ change. 

5.3 IMPLMEUIAIIOU 

State estimate based on real povver flows in lines and 
voltage magnitudes was obtained for IEEE 14 bus test system 
with gross error introduced in one of the real power flow 
data. The detection and identification of the presence and 
location of the had data were done by ,J(X) and r^j., tests 
respectively. The variances for all the measucements was 
taken as 0.01 p.u. and seven redundant measurements were 

made. 

Eetec t ion : 

J(X) test was applied to detect the presence of had 
data. The threshold value of J(X) was obtained from chi-square 
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distribution tables corresponding to 7 degrees of freedom 
and false alarm probability b = 0,005. Ibe test eould 
successfully detect the presence of gross measurement error 
beyond 15 The Pig. 5.1 shows the variation of 1(2) with 

magnitude of error. 

Identification : 

Identification of the bad data location was done 
for various cases both by r^j- and r^ tests. The r^^ teat 
always gave correct results but r^ test in few cases failed 
to identify correctly. A bad data having been detected 
and identified was removed from the input data, and the 
new J(X) was found to be well within limits of the chi-square 
value with reduced degrees of freedom and false alarm 
probability b = 0.005. 
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distrilDution tables corresponding to 7 degrees of freedom 
and false alarm probability b = 0.005. Ihe test could 
successfully detect the presence of gross measurement error 
beyond 15 W. The Pig. 5.1 shows the variation of J(2) with 
magnitude of error. 

Identification : 

Identification of the bad data location was done 
for various cases both by and tests. The r^^ teat 
always gave correct results but r^ test in few cases failed 
to Identify correctly, a bad data having been detected 
and identified was removed from the input data, and the 
new J(X) was found to be well within limits of the chi-square 
value with reduced degrees of freedom and false alarm 
probability b = 0,005. 
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OHAPISR 6 

33YRMIG state BSTBIAIIOR 


It IS often the case in closed loop control systems 
that the state variables which decide the optimal feedback 
control law are not directly measurable. Such missing 
state variables are to be reconstructed from the measurable 
output variables by observers. Pioneering work in the theory 
oi line.r time invariant deterministic and continuous time 
system observers was done by luenberger in 1964 [ 6 ] . Ihe 

theory has further been extended to nonlinear time tavariant, 
stochastic and discrete time systems later [55 , 57 , 41 ] 

6 . 1 APPLIOr-IIOF OP OBSSRYER THEORX 10 POWER 9 ISMS 
' Concept of such an observer theory to apply optimal 

control to synchronous machines has been drawing attention 
of many research workers. Significant improvement in the 
transient performance of power system can be achieved by 
applying optimal control law to synchronous machine exciters 
and governors. An observer can be built to asses the state 
of a synchronous machine and continuous or discrete feedback 
control signals can be generated from the knowledge of state 
vector. 

Optimal control of a linear dynamical system : 

A linear dynamical system is mathematically written 

as 
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X = A X + B u 

V - ' 

y = u X 

where x is n. state vector, u is r control vector and y is 
m measurement vector; A is nxn system matrix, B is nxr 
distribution matrix, 0 is mxn matrix which relates measure- 
ment vector y to the state vector x. 

The optimal control law in general is obtained for 
a systom of above kind by choosing u such that the perfor- 
mance index 

J='^ / (y Qy+u^R u)dt 
0 

is minimized, where Q and R are the weighting matrices. 
Solution u IS obtained by minimum principle 

u (t) = -R ' B K x(t) 

where K is symnetric positive definite solution of the 
matrix Riccati equation 

KA + A^(K-K B R“^ B^ K + 0 Q G) = 0 

The state vector x(t) can be obtiined from an observer. 

In [lO] Arumugam has shown improvement in the 
transient behavi'-'ur of a salient pole machine connected to 
infinite bus by applying such an optimal control law to 
excitersand governors using a linear observer. The dynamical 
model of the synchronous machine was linearized to the 
standard form of (6.1). 
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Such a linearized model approach may not always "be 
permissible because the observer of the linearized model 
IS very sensitive to the initial starting value of the 
unknown state vector x. Also as there is a growing trend 
towards the use of discrete supplementary excitation control 
during transient period and use of continuous or discrete 
control in dynamic state of the synchronous machine, the 
transient state needs “to be accurately sLasessod which a 
linearized model may fail to do successfully. 

An extended linear observer has been proposed in [ 23 ] 
to assess the state of a synchronous machine more accurately, 
'•^he observer has been shown to be less sensitive to the 
starting value of the state vector compared to a linear 
observer given by Arumugam in [10] . The observation 
vector y is obtained by a half cycle mean value detecting 
circuit given in [24] . This circuit holds the average value 
over half power frequency cycle of the input quantity at its 
output terminals, 

6.2 EXTENDED LINEAR OBSERVER 

Nonlinear power system equations are. written as 

x(t) = f(x(t) ) (6 » 2 ) 

where x is n state vector andf(») is a nonlinear vector 
valued function. The m size observation vector y is obtained 
sequentially at every power frequency cycle by the mean value 
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detecting circuit of [24]. Ihe nonli near relation "between 
y and X IS 

(6.3) 

where K = 1 , 2 ,.... corresponds to period of the power 
cycle T, The estimated value of state vector at 

interval according to linear observer theory is given 

by 

^ ^K/xC-l “ ^K/K~1^ (6*4) 


iC/lC-l predicted value of x^- from the previous 


where x 

estimate at (K-1)th interval, -| the predicted obser 

vation value of y-^ from ^^i/YL 1 ^ suitably chosen 

filter matrix of nxm size, is the measurement vector at 
kth interval, 

^K/‘C-1 obtained by Taylor series expansion upto 
second order term as 




T 


“ ^ 4-1 ^ * K -1 ^-1 


x=x, 


■K-1 


m2 . 

+ -2- f I _ fj (6.5) 

^=^-1 ^-^-1 


'The value of obtained from (6.5) at Xg._^ is written 
as ^K/K 1 value of y-^ obtained by substituting 

%/iC-1 (6.3). The filter matrix % which has the 

properties desired for convergence of eqn. (6.5) is put as 

% = (%%; + 


( 6 . 6 ) 
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where 


and 






[ 0^] = 0.01 [Hj ] 

X ^ JL 

[ Cg] =0 (i ^ 3) 

X * 3 


(6.7) 


( 6 *S) 


6.5 im'LEMEFMION 

The above extended linear observer was tested for 
the cylindrical rotor machine connected to infinite bus given 
[ 23] • actual transient response was obtained by 

fourth order H-K method. The results of the muthod were 
given as input measured quantities y^ to the extended linesr 

observer. The response of the obsoryer is shown in 
Eig. 6.1. e 

The observer described m this chapter applies to 
uncontrolled case. It remains to be demonstrated how this 
observer can be modified for a controlled case. 
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EURTHER RESEARCH TOPIOS IN STATE ESTIMATION 

There are two ma^or class of problems associated with 
state estimation. Eirst results from the difficulties 
encountered in carrying out the control applications of' 
stati c and dynamic estimation successfully. Second problem 
IS related to improvement of the static state estimation 
accuracy . 

7.1 TOPICS related to STATE ESTIMATION APPLICATIONS 

Many potential applications of state estimation to 
improve static and dynamic performance of power system 
discussed in Chapter 1 may become offset because of the 
following practical difficulties: 

(1) In an interconnected system the information about the 
external system is lacking in many respects (network topology, 
line-parameters, voltage conditions etc.). 

(2) Machine models for transient state are inaccurate 
in the sense that system nonlinearities like saturation etc. 
are not taken into account by the models. Model parameters 
change with changing load generation pattern. 

The first difficulty offsets the functions of steady 
state security assessment, contingency evaluation and stability 
investigation for a simulated critical continfc,ency in own 
system, whereas the second ''difficulty impairs the performance of 
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the optimal controllers used for a synohronous machine. 

Each of these difficulties form a different class’ of 
problems. 

External system equivalents for steady state security 
assessmen t 

Ihe external system equivalencing for steady state 
security assessment is basically a problem of finding 
admittance parameters of the reduced equivalent external 
system. The external system is reduced to managable size 
and augmented to the own system (O.S.) such that the 
"buffer zone" is a true representation of external system 
(E.S.) for a particular load generation pattern and status 
of the external network. Very few research workers have 
looked into this problem. Debs is notable among them [16,22]. 
Debs has recently proposed what seems to be a satisfactory 
solution to this kind of equivalencing problem [ 3 I]. 

Dynamic equivalents 

It is desired to predict the effect of a critical 
contingency in own system on stability of the overall system. 
The external system is represented by reduced order dynamic 
system by lumping the coherent group of generators together 
in the external system. The coherency analysis and dynamic 
equivalence is a topic of much research today . 
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or realization, adaptive control : 

The second category of difficulties have to be 
overcome by applying adaptive control techniques. The 
parameters of model under consideration are identified or 
realized and the mod< 1 adapted continuously. Application 
of such an identification problem can be found in [ 19 ]- 

7.2 TOPICS RELATED TO IMPROYMMT II THE EEEICIEYCY OE 
STATIC STA.TE ESTIMTORS 

Efforts have been made to improve the static staf’e 

estimators efficiency in terms of accuracy of results, 

(1 ) Choice of meter locations in the system: 
ce of meter locations in the system : 

The question of selecting the best choice of meter 
location to improve the quality of results is an important 
one ■which has not been answered satisfactorily so far. 

However, in [ 20] some attempt has been made to answer the 
question of deciding the best measurement to add to an 
existing measurement set. It follows from linear estimation ^ 
theory of statistics that the best choice of the matrix H 
relating observation vector ^ with state vector X as Z = H X wil 
will be to have columns of H orthogonal to each oth^'r under ^ 
the condition that = constant, where is the ith t 

column of H (]40] » Selecting H is not in our hands for 
power system problems but one could select the best measurefient , 
to add to a set of existing measurements by checking ortho- j 
gonality of colu ms of H for each of the choices of measurement^ 
at disposal. It is proved in [ 20] that most orthogonal vector j 
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to the coluans of a matrix H corresponds to eigenvector 
associated with the smallest eigenvalue of the matrix 
Use IS made of the above theorem in [20] to decide best 
measurement to add to a set of measurements but the [2*14] 
computations involved are prohibitively large. A considerable 
amount of work needs to be done in this area. 

K alman filtering ; 

Recursive xCalman filters have been applied to 
state estimation by modelling the power system as 

’Si] ^ 

where w^^ is a stochastic process of known statistics. This 
needs a statistical load modelling of the power system loads 
at different intervals of time. Kalman filtering solution of 
state esti nation is more realistici in the sense that the 
power system is represented by dynauiical equations. Accurate 
knowledge of statistics of the stochastic process is a 
pre-requisite for Kalman filter solution. 
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OHAPm 3 

EVALUA'IIQTT 0? DIPF31TIEITT ALGORITHMS klVD GQT^LUSION 

The static state estimation algorithms discussed in 
Chapters 2, 5 and 4 are compared in this chapter in terms 
of computational time, storage requirements and their merits 
and limitations. Each of these methods is evaluated below. 

8.1 I^TEElETOEtTT E^UATIOTT SOLUTIOTT METHOD 

It is seen from the Table 3.1 that this method 
requires minimum time and storage (In fact the storage 
requirements can further be reduced as there were many 
variable names used fou the sake of clarity of the programme). 
Basically this method gives an approximate solution to state 
estiaation and accuracy of results is largely dependent upon 
system configuration. In spite of the Urn tation this method 
has, there is a distinct advantage of this method over others 
that it can be implemented on a small computer. It does not 
involve matrix storage and trigonometric functions. Hence 
memory size and system software and hardware for the 
computer 3.r6 mininiuiii# 

This method assumes that line flows at both ends 
are available. However, if one end line flow is missing 
estimation should he carried out with only one end line 
flow instead of putting the missing end line flow equal to 

Alternatively if th e bus injection 


the other end line flow. 
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measurement at the missing end bus is available the missing 
line flow can be reconstructed and used for estimation^ Sus 
infection if available could also be used to check which end 
line flow measurement is to be rejected when the two end 
line flows are found to be significantly different from 
each other and one of them is thought to be incorrect. 

An improvement in the results could possibly be 
obtained by taking average of the solutions obtained by 
tracing alternative paths. It is of course not possible 
to trace all possible paths and to obtain all possible sets 
of solutions but a properly selected group of paths will 
give good enough results. 

B.2 ’LO' AIGORITHM 

This method is seen to have taken next lowe:^rmemory 
storage and computer time. The only approximation involved 
in derivation of this algorithm is that the weights to the 
P and <4 measurements at one end of line flow are assumed to 
be equal. 

The mam advantage of this algorithm is that it is 
very easy to build the gain matrix and to modify it for 
changes in system configuration. The gain matrix is constant 
and sparse in nature. When e ual weights are given all 
measurements, i.e. W is taken as identity the off-diagonal 
entries of the gain matrix become 0, 1 or -1 . Individual 
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memory bits can be utilized to store these entries instead 
of storing them m a single word. Such storing and retrieving 
of individual bits can be done by AUD and OR functions of 
PO'HRan 'Vvliich. are avah-lable in most of the modern coia,r,ut ers . 

When these techniques for efficient storage are applied to 
’ LO’ method, it becomes more powerful, iloreover the time 
talcen to converge given in fable 8.1 for this method is 
rather misleading, it will take lesser time in on line 
mode of operation. The gain matrix does not change unless 
the system configuration or meter location itself changes.' 

Thus the inversion or triangular factorization of the gain 
matrix which is the most time consuming part of the algorithm 
IS obviated in real time mode of operation. 

It IS assumed in this method that the reference bus 
voltage IS correctly known. The error in it affects the 
results. 

S.3 ’PV’ METJOD 

Thi s method can be seen from Table 8.1 to have 
moderate memory and time requirements. This method is not 
different from a lumped method. The objective to test 
this method was to demonstrate the feasibility of solution 
from real line flow and bus voltage measurement data. This 
method is particularly suited for Indian conditions as the 
data it takes is readily available and the additional cost 
of equipment is not required. 
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S . 4 ' FZDS ' AND ' LUMPED WIS' ALGORITHM 

Past decoupling approach basically is used to reduce 
the memory and convergence time of lumped WLS method. Prom 
the Table 8.1, it can been that PDS does take lesser time 
than a lumped ¥LS method but no such improvement in total 
memory storage is apparent. This is because of the fact 
that the sample system on which these methods have been 
implemented is very small in size. The data storage in lumped 
’JLS is api/roximately proportional to n for an n bus system to 
n^/2 for PDS method. 

Unlike 'LO' method, PDS and lumped WLS methods can 
take any number of measurement and a very high figure of 
redundancy can be obtained with them. 
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Table 8.1 


Comparison of different methods. 


’LO’ Lumped WLL EDS P7 Independent 

method method method method equation 

solution 


;a 

^ ' oraf^e 

1.7 K 

5 K 

Programme 

storage 

1.3 K 

3.9 K 

To tal 
storage 

3.5 K 

3.9 K 

Time talcen 

7.3 sec. 

133 sec. 

Total no. of 
iterations 

5 

4 


7.7 K 

3.6 K 

0.296 K 

4.2 K 

2 K 

0.924 K 

9.9 K 

5.6 K 

1.22 K 

99.66 

sec. 36,65 

sec. 2,466 sec 

12 

4 

Max. 3 
iterations 
per line 
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LEAST SQUARES ESTBIATIOR: 

Consider a static system in which an observable 
variable vector z of size m is connected with n state vector x* 
z = H X (A.1) 

where H is a known mxn matrix. It is desired to solve for 
X from given z and H. Assume that m^n. If all the observa- 
tions were to be correct, eqn.(A.I) would be a set of consistent 
equations, i.e. any n equations out of the m would have given 
the same solution to x. However if(A.1) is not a set of 
consistant equations, the vector x is chosen which minimizes 
the performance index 


J = (z- H x)^ (z- H x) 


T 


or dJ/dx= (-H) 2(z- H x) = 0 


X 


m -1 1 

= (H H) H z 


(A.2) 


(A,3) 


WEIGHTED LEAST SQUARES ESTIJ-IATIOF: 

Consider a case where the observations z are known to 

be contaminated by a random noise 


z = H X + h 


(A. 4) 


Estimation of state x is done Dy the minimisation of weighted 
aum of squares of errors by giving weights to each measurements 


of z, as follows 



%z 


where ¥ is the inverse of error covariance of the noise 
It becomes a diagonal matrix when the noise vector is 
uncorrelated. 

The stimate of x can be obtained similar to(A.3) 
as follows 

X =(H^¥ H)“'' ¥(z-Hx) (A.6) 
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IEEE 14 BUS TEST SYSEEM 


Ihe 14 Bus test system shown in Pig. 

TaBle B.1: Impedance and line 
charging data# 


Line 

Desi gnation 


Resistance 

p.u. 


Reactance Line charging 


p.u. 


p.u. 


1 

1 

2 

2 

2 

3 

4 
4 

4 

5 

6 
6 
6 
7 
7 
9 
9 

10 

12 

13 


- 2 

- 5 

- 3 

- 4 

- 5 

- 4 

- 5 

- 7 

- 9 

- 6 
- 11* 
- 12 

- 13 

8 

- 9 

- 10 

- 14 

- 11 

- 13 

- 14 


0.01933 

0.05405 

0.04699 

0.05811 

0.05695 

0.06701 

0.01335 

0,0 

0.0 

0.0 

0.09498 

0.12291 

0.06615 

0.0 

0.0 

0.03181 

0.12711 

0.08205 

0.22092 

0.17093 


0.05917 

0.22304 

0.19797 

0.17632 

0.17388 

0.17102 

0.04211 

0.20912 

0.55618 

0.25202 

0.19390 

0.25531 

0.13027 

0.17615 

0.11001 

0.08450 

0.27038 

0.19207 

0.19933 

0.34802 


Impedance and line charging susceptanoe 
on a 100,000 kva Base, 
line Charging one-halfof total charging 


0.0261 
0.0246 
0.0219 
0.0187 
0.0170 
0.0173 
0.0064 
0.0 
0.0 
0.0 
0,0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

is per unit 
of line. 
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Table B.2; Generation and load ScheduleB, 


Bus 

Number 

Generation 

load 

MW 

M7AR 

MW 

MVAR 

1* 

0.0 

0.0 

0.0 

0.0 

2 

40.0 

0.0 

21.7 

12.7 

3 

0.0 

0.0 

94.2 

19.0 

4 

0.0 

0.0 

47.8 

-3.9 

5 

0.0 

0,0 

7.6 

1.6 

6 

0.0 

0.0 

11.2 

7.5 

rj 

0.0 

0.0 

0.0 

0.0 

S 

0.0 

0.0 

0.0 

0.0 

9 

0.0 

0.0 

29.5 

16.6 

10 

0.0 

0.0 

9.0 

5.8 

11 

0.0 

0.0 

3.5 

1.8 

12 

0.0 

0.0 

6.1 

1.6 

13 

0.0 

0.0 

13.5 

5.8 

14 

O.o 

0.0 

14.9 

5.0 

* indicates slack bus 

» 




Table B.3: Transformer lata. 


Transformer 

designation 

4-7 

4- 9 

5- 6 


Tap Setting 


0.978 

0.969 

0.932 



Table B.4: Static Capacitor Data 


Bus Nuaber 

Susceptance p.u. 

9 

0.19 


Table B.5: Load Plow Solution 


Ib?om 


Line Plow 
M¥ M7AR 


156.8 

75.6 
-152.5 

73.2 

56.2 

41.5 

-70.9 

-23.3 

-54.5 

23.7 
-61.2 

28.0 

16.1 

-72.8 

-40.6 

61 .7 
44.0 

-44.0 

7.3 

7.3 

17.7 

-28.0 

0.0 

^. 8.1 

0.0 

-16.1 

-28.1 

5.2 

9.5 
-5.2 
-3.7 
-7.4 

3.7 

-7.7 

1.6 

-17.5 

- 1.6 

5.6 

-9.3 

-5.6 


-20.4 

3.5 

27.6 

3.6 
-2.3 

0.7 

1.6 
2.8 

3.4 
-5.4 

15.6 
-9.4 
-0.3 

2.6 

- 1.6 

-15.3 

12.8 

-8.4 

3.5 

2.5 
7.2 

11.1 

-16.9 

6.8 

17.4 

1.6 

-5.0 

4.4 

3.7 
-4.3 
- 1.6 
-3.4 

1.6 

-2.4 

0.7 

- 6.8 

-0.7 

1.7 
-3.4 
- 1;6 



Table B.5 (Continued) 


Bus 

Voltage 

Generation 

Load 

Magnitude 

Angle 

m 

MVAR 

m 

MVAR 

1 

1.06 

0.0 

232.4 

-16.2 

0.0 

0.0 

2 

1.045 

-4.98 

40.0 

42.4 

o 

♦ 

CM 

12.70 

5 

1.01 

- 12.72 

0.0 

33.4 

94.2 

19.00 

4 

1.019 

-10.33 

0.0 

0.0 

47. S 

-3.9 

5 

1.02 

-8.78 

0.0 

0.0 

7.6 

1.6 

6 

1.07 

-14.22 

0.0 

12,2 

11.20 

7.50 

7 

1.062 

-13.37 

0.0 

0.0 

0.0 

0.0 

8 

1.09 

-13.36 

0.0 

17.4 

0.0 

0.0 

9 

1.056 

-14.94 

0.0 

0.0 

29.50 

16.60 

1 0 

1.051 

-15.10 

0.0 

0.0 

9.0 

6.8 

1 1 

1.057 

-14.79 

0.0 

0.0 

3.50 

1.80 

12 

1.055 

-15.07 

0.0 

0.0 

6.10 

1.6 

13 

1.05 

-15.16 

0.0 

0.0 

13.50 

5.8 

14 

1.036 

-16.04 

0.0 

0.0 

14.9 

5.0 


miA INPUT TO VARIOUS STATE ESTIMATION ALGORITHMS: 

1. Lumped ¥LS and EES Method 

Bus injections at all buses, line flows in one 
direction of all lines and voltage magnitudes at all buses 
were used as input data for both lumped WLS and EDS methods 
Eollowing is the description of lines considered. 
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LINE FLOWS 

From 

To 

From 

To 

Fro 

To 

From 

To 

1 

2 

3 

4 

6 

11 

9 

10 

1 

5 

4 

5 

6 

12 

9 

14 

2 

5 

4 

7 

6 

13 

10 

11 

2 

4 

4 

9 

7 

8 

12 

13 

2 

5 

5 

7 

7 

9 

13 

14 


2. PV Method 

The real power flow in the lines (one end only) 
considf rod in the PDS and lumped WIS method and voltage 
magnitudes at all buses were given as input data to PV 
algorithm. 


3. LQ Method 

Line flows in both the directions of the following 
lines was the input to this algorithm. 


lines 

From 

To 

• 

From 

To 

From 

To 

From 

To 

1 

2 

4 

7 

7 

9 

13 

14 

1 

5 

5 

6 

9 

10 

4 

5 

2 

3 

6 

12 

10 

11 

2 

4 

3 

4 

7 

3 

12 

13 

9 

14 
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4. Independent Equation Solution ; 

Ihe line flows in "both directions of the following 
lines was the input data. 


From 

1 

2 

2 

2 

4 

4 

9 

9 

5 

6 

6 

6 

7 

To 

2 

3 

4 

5 

7 

9 

10 

14 

6 

12 

13 

11 

8 
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APPENDIX G 

DESGRIITION OF PHE SYNCHRONOUS mOHINE USED FOR TRANSIENT 

SOIATE ESTIMTION PURPOSE 


Phe following assumptions were made for the synchronous 
machine used for transient state estimation studies agcby the 
extended linear observer. 

1 . Cylindrical rotor with no damper windings. 

2. Pwo fields with same structure placed on the d-and q-axis. 

3. No saturation 

4. Transformer voltage, armature resistance and line 


resistance neglected. 



Synchronous machine connected to infinite bus 
Nom enclature 


T - field open circuit time constant 

do 

X - synchronous areactance 

s 

X^ - line reactance 

T - field short circuit time constant 

dl 

e & e~ - direct and quadrature axis component of the 
f d f <1 

field excitation voltage 


w 

P 

in 

H 


- angular speed rad/sec 

- mechanical power input 

- inertia time constant (sec.) 
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moment of inertia (= 2H/w) 
damping coefficient. 


Model of the Synchronous Machine 
State vector x consists of 
X = (x^ , X2» -siy = (e^, 6 , 6 

She model x = f(x) is -written as 

= -A^x^ + A2 X^cos x^ + A^sin x^ - 

Xg = “■^1^2 ~ ^3 ^3 + -^5 

X3 = X4 

x^ = A^x^ + Ag(x2 sin x^ - x-jCos x^) + Ag 
where 

^o = "do -^"dV ^1 = ' 

^2 = "d1 ®s/\ ’ S " ’ 


( 0 . 1 ) 


( 0 . 2 ) 


(0.3) 




A = e. /A ; 
5 fq 0 


A, = P. /M ; 
6 in 


-D/M ; 


8 X. 


She output vector y consists of terminal voltage y^, output 
power yg speed deviation y^. She relation y = h(x) is 


written as: 

t 2 2.^ 

y^ =(x^ + Xg) 


^2 


(xg sin X3 - x-| cos X3) 


(0.4) 



Systrm constants 


Xg = 1, x^ = 0.6, = 3 sec., = 1 sec. 

H =5 sec., D = 0.005, w = 2n(60) rad/sec., 


®fd“ ^in “ 0.76923. 


Initial values of the state vector 


= -0.48595, = 0.86149, = 0.507635 
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